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Abstract. We present sub-arcsecond images of AB Auriga obtained with the IRAM Plateau de Bure interferometer in the 
isotopologues of CO, and in continuum at 3 and 1.3 mm. These observations allow us to trace the structure of the circumstellar 
material of AB Aur in regions where optical and IR imaging is impossible because of the emission from the star. These images 
reveal that the environment of AB Aur is widely different from the proto-planetary disks that surround T Tauri stars like DM Tau 
and LkCal5 or HAeBe stars like MWC 480, in several aspects. Instead of being centrally peaked, the continuum emission is 
dominated by a bright, asymmetric (spiral-like) feature at about 140 AU from the central star. Little emission is associated with 
the star itself. The molecular emission shows that AB Aur is surrounded by a very extended flattened structure ("disk"), which 
is rotating around the star. Bright molecular emission is also found towards the continuum "spiral". The large scale molecular 
structure suggests the AB Aur disk is inclined between 23 and 43 degrees, but the strong asymmetry of the continuum and 
molecular emission prevents an accurate determination of the inclination of the inner parts. Analysis of the emission in terms of 
a Keplerian disk provides a reasonable fit to the data, but fails to give a consistent picture because the inclinations determined 
from 12 CO J=2->1, 13 CO J=2^1, 13 CO J=1^0 and C ls O J=1^0 do not agree. The mass predicted for the central star in 
such Keplerian models is in the range 0.9 to 1.2 M Q , much smaller than the expected 2.2 M Q from the spectral type of AB Aur. 
Better and more consistent fits to the 13 CO J=2 — > 1 , 13 CO J=l— >0 data are obtained by relaxing the Keplerian hypothesis. We 
find significant non-Keplerian motion, with a best fit exponent for the rotation velocity law of 0.41 ± 0.01, but no evidence for 
radial motions. The disk has an inner hole about 70 AU in radius. The disk is warm and shows no evidence of depletion of CO. 
The dust properties suggest the dust is less evolved than in typical T Tauri disks. Both the spiral-like feature and the departure 
from purely Keplerian motions indicates the AB Aur disk is not in quasi-equilibrium. Disk self-gravity is insufficient to create 
the perturbation. This behavior may be related either to an early phase of star formation in which the Keplerian regime is not 
yet fully established and/or to a disturbance of yet unknown origin. An alternate, but unproven, possibility is that of a low mass 
companion located about 40 AU from AB Aur. 
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1. Introduction 

The existence of disks around Herbig Ae/Be stars is now 
widely accepted but the physical properties of these disks are 
poorly known. Contrary to TTauri disks, there are only a few 
examples of large Keplerian disks a round Herbig Ae stars: th e 
A4 star MWC480 JMannings et alJI 19971 Isiinon et alJl2000l) . 
and the AO star HP 34282 JPietu et all2003T) ~ 

AB Auriga is one of the nearest, brightest and best stu died 
Herb ig Ae stars. It has a spectral type A0-A1 ( iHernandez et alJ 
2004) and is located at a distance D = 144+^ pc following 
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Hipparcos measurements Jvan den Ancker et alJ fl9 98). NIR 
imaging with the HST-STIs lGradv et alJll999T reveal a large 
envelope surrounding the star and scattering the stellar light. 
This flattened reflection nebulae is seen close to pole-on, up 
to r ~ 1300 AU from the star. New NIR observations, per- 
formed with the Subaru telescope us ing the Coronographi c 
Imager and Adaptive Optics systems (Fukagawa etalJ l2004). 
show that the circumstellar matter presents a spiral structure. 
AB Auriga has been also observed in MID-IR spe ctroscopy 
with ISO JMeeus etalJl200ll iBouwman etaD l2000l) . These 
data constrain the dust content several tens of AUs from the 
star. The modelling of the SED in this frequency range als o 
shows that the star belongs to the Group I jMeeusetalJ200ll) . 
In the classification of Herbig Ae circumstellar matter, stars of 
Group I are surrounded by a flaring disk responsible for the 
bum p in the M ID-IR. The disk was imaged at 1 1 and 18yum by 
IChen & Jural §003), using the Keck I telescope. They found 
the disk is resolved at 18.7/vm with an approximate diameter 
of ~ 1.2" or 170 AU. Technical progresses achieved on op- 
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tical interferome ters allow the first quantitative studies of the 
very inner disks. M onmer &Millan-Gab"etlll2002h have shown 
that in such objects, the very inner dust disk is truncated by 
sublimation of the dust at temperature around ~ 1500 K. For 
AB Auriga which is a star of r e (j ^ 10000K, this happens at 
r,„ ~ 0.2 - 0.3 AU. Finally, the surround ing material was also 
observed at millimeter wavelengths by Mannings & Sargent 
il 19971) with OVRO, but the derived inclination 75°) is much 
larger than that found from all other (optical or IR) tracers 
O 20 - 30°). 

Following our successfu l study of DM Tau in CO isotopo- 
logues (Darto iset alJl2003h . we decided to perform a similar 
study of the AB Auriga environment. Therefore, we observed 
AB Auriga with the IRAM interferometer in I2 CO J=2— >1, 
13 CO J=1^0 and J=2^1, and C I8 J=1^0. We report here 
the results of this multi-line, CO isotope analysis. Our goals 
are double: we want to 1) characterize the physical structure of 
the circumstellar matter (vertical temperature gradient, density, 
kinematics) and 2) compare quantitatively the large scale mm 
properties with those of TTauri disks. The observations and the 
results are described in Sections 2 & 3, respectively. Section 4 
presents the best model and we discuss in Section 5, the physi- 
cal implications. 

2. Observations and Results 

2.1. PdBldata 

The 12 CO observations used 5 antennas and were carried out 
in winter 2001/2002 in D and C2 configurations. Baselines up 
to 170m provided 2.00 x 1.60" resolution at PA 92° for the 
1.3 mm continuum data. We observed simultaneously at 89.2 
GHz (HCO + J=l— >0, which will be discussed in a forthcom- 
ing paper) and 230.5 GHz ( 12 CO J=2->1). At 1.3 mm, the tun- 
ing was double-side-band (DSB) while at 3.4 mm, the tuning 
was purely single-side band (LSB). The backend was a cor- 
relator with one band of 10 MHz (spectral resolution 0.23 
km.s 1 ) centered on the HCO + J=l— >0 line, one band of 20 
MHz (0.18 km.s~' resolution) centered on the 12 CO J=2->1 
line, and 2 bands of 160 MHz for the 1.3 mm and 3.4 mm 
continuum, respectively. The phase and flux calibrators were 
0415+379 and 0528+134. The rms phase noise was 8° to 25° 
and 15° to 50° at 3.4 mm and 1 .3 mm, respectively, which intro- 
duced position errors of < 0.1", and a seeing better than 0.3". 
The observation time was shared with other sources: MWC480 
and LkCal5 or CQ Tau and MWC758, AB Aur using only 
a small fraction of the transit time. The total on source inte- 
gration time is ~ 10 hours. As a consequence, the amplitude 
and phase calibration of these sources is homogenous and the 
flux density of all sources were estimated simultaneously. This 
allows us to make accurate comparisons of the spectral in- 
dex of the continuum emission. In all cases, the flux density 
scale was referred to MWC349, for which we used a flux of 
S(v) = 0.95Jy(v/100 GHz) 06 . 

The C 18 OJ=1^0, 13 COJ=1^0 and J=2^1 observations 
were performed between 2001 and 2004. Configurations D, 
C2, B 1 (5 antennas) and A (6 antennas) were used, and pro- 
vided baselines up to 400 m. This provides angular resolu- 



tion of 0.80 x 0.55" at PA 39° at 220 GHz, with a residual 
seeing < 0.2" and position errors < 0.05" after calibration. 
The correlator provided a spectral resolution of 0.09 km.s 1 
for 13 COJ=2^1 and O.lSkm.s" 1 for the 13 COJ=1^0 and 
C ls O J=l— »0 lines. The same phase and flux calibrators were 
used. The total integration time is around 30 hours, half of 
which was spent on the long baseline (A configuration) obser- 
vations. 

We used the GILDAS software package to reduce the data. 
Images are presented at different angular resolutions, obtained 
by applying a taper and re-weighting the data when needed. 
However, the natural weights were used in the disk modelling. 

2.2. 30-m data 

In 2000, we also obtained a single-dish spectrum of 
13 CO J=2— >1 at the 30-m telescope, which when compared to 
the PdB measurements indicates at most 20 % of the flux was 
missed by the interferometer. Accordingly, it is unlikely that the 
13 CO measurements are significantly affected by missing flux, 
specially within the inner 10 to 20". However, the 12 CO J=2— > 1 
line is much brighter and most likely hampered by structures 
within a more extended region. 

3. Results 

Figure [2 is a montage displaying high resolutions images of 
the continuum emission at 2.8 mm and 1.4 mm, and of the 
integrated line emission of 12 COJ=2— >1 13 COJ=2— >1 and 
I3 COJ=l— >0 transitions. The 1.4 mm emission was obtained 
by merging the 220 and 230 GHz data in together to improve 
signal to noise; the flux was scaled to a reference frequency of 
220 GHz assuming a spectral index of 3.0 for the emission. The 
(lower resolution) emission at 1 10 GHz is also presented. The 
total continuum flux at 1.4 mm is 85 ± 5 mJy, and at 2.8 mm, 
11 ±2 mJy. 

Figure |2 presents the channels maps of the 12 COJ=2— >1, 
13 COJ=2->l and 13 COJ=1->0 transitions. The 13 COJ=2^1 
data has been smoothed to about 1" resolution to help the com- 
parison with the J=l— >0 transition. 

In Figure [3] the co ntinuum emi s sion is super-imposed on 
the STIS image from iGradv et af] dl999l) . No position ad- 
justment has been made in doing so: the absolute astromet- 
ric accuracy of the PdB measurement is better than ^ 0.05". 
The best fit position of the centroid of the 13 CO emission is 
a = 04:55:45.843, 6 = 30:33:04.21 only 0.03" from the extrap- 
olated Hipparcos position. The HST STIS image (which was 
centered on the Hipparcos coordinates for 1999) was shifted 
southward by 0.08" to correct for the proper motion of AB 
Aur (taken from the Hipparcos catalog) between the epochs of 
the two observations (1999 and 2003). The emission is clearly 
concentrated in an arc-like structure about 1" from the central 
star. The flux density measured at 1.4 mm in the direction of 
the central star is ^ 5 ± 1 mJy, but this should be considered 
as an upper limit because of contamination from the arc-like 
structure due to the limited angular resolution. 

In Figure n the line emission was integrated between 3.8 
and 8.2 km.s . This interval contains most of the line pro- 
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Fig. 1. Line emission (integrated over the [3.8,8.2] km.s -1 velocity range, and continuum emission towards AB Aur. The synthe- 
sized beam is indicated in each panel. The cross indicates the direction of the major and minor axis of the AB Aur disk, as derived 
from the analysis of the line data, but using an inclination of 30° for the aspect ratio of the cross. Top row, left: Integrated intensity 
of the 12 CO J=2— >1 line. The angular resolution is 2.0 x 1.6" at PA 92° and the contour spacing is 0.85 Jy/beamkm.s -1 or 12 cr. 
Middle: Integrated intensity of the 13 CO J=2— >1 line. The angular resolution is 0.72 x 0.63" at PA 29° and the contour spacing 
is 85 mJy/beamkm.s 1 or 2.2 cr. Right: Integrated intensity of the I3 CO J=l— >0 line. The angular resolution is 1.37 x 1.05" at 
PA 25° and the contour spacing is 34 mJy/beamkm.s -1 or 2.3 cr. Bottom row, left: continuum emission at 2.8 mm. The angular 
resolution is 1.37 x 1.07" at PA 27° and the contour spacing is 0.4 mJy/beam, 28 mK, or 2.0 cr. Middle: continuum emission 
at 1.4 mm. The angular resolution is 0.85 x 0.59" at PA 18° and the contour spacing is 1 mJy/beam, 51 mK, or 2.7 cr. Right: 
Integrated intensity of the C ls O J=l— >0 line. The angular resolution is 1.37 x 1.05" at PA 25° and the contour spacing is 25 
mJy/beamkm.s 1 or 2.1 cr. Coordinates are J2000.0 Right Ascension and Declination. 



file, but weak emission in the line wings is still detectable 2 
km.s away (specially in I2 CO). The integrated line emission 
maps reveal a clear inner hole in the distribution of the emis- 
sion from the CO isotopes, as well as enhanced emission at the 
location of the continuum ring. Note in particular that, while 
the 2.8 mm continuum hardly shows the central depression, the 
13 CO and C I8 OJ=l— >0 line emissions obtained at the same 
angular resolution clearly show the ring. 

3.1. Analysis Method 

To better quantify the properties of the AB Aur circumstellar 
materia l, we used an appro ach similar to that applied to DM 
Tau bv iDartois et alJ y003). The analysis is based on a new 
radiative transfer code, with non-LTE capabiliti es, whi ch will 
be presented in a forthcoming paper ( Pietu et alJ l2005h . Apart 
from discretization details, the LTE part (which was used ex- 
clusively here) is however equivalent to the code described 



in lDutrev et all d 1994 . As in lGuillotea u & Dutrev ( 1998), the 
comparison between models and data and the minimization is 
done inside the UV plane, using natural weighting. The physi- 
cal model is that of a rotating, flared disk in hydrostatic equilib- 
rium, with power law distributions for the density, temperature, 
rotation velocity and scale height as function of radius. Since 
the AB Aur disk is relatively extended, we have applied pri- 
mary beam attenuation to the model images before computing 
the model visibilities. We have also improved the minimization 
technique by searching the minimum of the x 1 with the help of 
a Levenberg-Marquard method, instead of a simple grid search 
previously. The two methods were intensively compared, the 
new one improving the convergence speed by a factor ~ 20-30 
on typical data, and providing more accurate estimate of the 
error bars since it takes into account the correlation between 
parameters. Table ^ summarizes the parameters found for the 
best models for the various lines and the continuum. 
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Fig. 2. Channel maps of the CO isotopologues in AB Aur. The LSR velocity, in km.s" 1 , is indicated in the upper left corner of 
each panel. Coordinates are in arsec from the star position. All data have been smoothed to the same spectral resolution, 0.25 
km.s . Top: 13 CO J=l— >0 line, spatial resolution 1.37 x 1.05" at PA 28°, contour spacing 30 mJy/beam, or 2.1 K, or 2.2 cr. 
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Fig. 3. 1.4 mm conti nuum data (in con tours) superimposed on 
the HST image from lGradv et alJ lll999l) . in grey scale. The an- 
gular resolution is 0.85 x 0.59" at PA 18°. The contour spacing 
is 0.9 mJy/beam, corresponding to 45 mK, or 2.4 <x, with the 
zero contour omitted. The cross is as in Fig. 1 



3.2. Continuum Emission: a spiral pattern? 

The continuum emission from AB Aur is quite peculiar. While 
most other T Tauri stars and HAeBe stars observed at high an- 
gular resolution show a centrally peaked emission, the emission 
from AB Aur is reminiscent from that of GG Tau, except for the 
angular scale which is a factor of 2 smaller here. The emission 
peaks accurately coinc ide with br ight features observed in the 
STIS image of lGradv et alJ dl999h . but our observations allows 
to trace the regions which are occulted by coronagraphic mask 
or by the telescope spider in the HST observations. 

No reliable estimate of the inclination of the structure can 
be obtained from the 1 .4 mm image only, because of the highly 
asymmetric nature of the emission. We have deprojected the 
emission assuming at position angle -31°, as derived from our 
analysis of the 13 CO data, and an inclination of 23°. These 
values are also in ag reement with the numbers quoted by 
iFukagawa et alJ J20 04) from the scattered light in H-band. The 
result is presented in Fig@] the asymmetry seen in continuum 
may follow i nwards the spira l -like f eatures detected in scat- 
tered light by Fukagaw a et alJ d2004l) . By extension, we shall 
refer to this structure as "the spiral". 

Although the spiral near 160 AU is the most striking fea- 
ture in the continuum, weak extended emission is also found 
at larger distances from the star. When interpreted in terms of 
a truncated disk, the best fit requires a disk with an inner ra- 
dius of =i 110 AU, an outer radius 350 AU, and a power 
law distribution for the surface density of p 2.3 (see Table 

This outer r adius is with i n 2cr of the 280 AU ring detected 
at 20.5 fim bv IPantin et alJ ll2005h . Negligible emission arise 
from beyond 350 AU, as shown by the lower resolution, higher 
brightness sensitivity images displayed in Fig|5] A larger inner 




Fig. 4. Deprojected image (from an inclination of 23°) of the 
continuum emission of AB Aur at 1 .4 mm. The black curves 
indi cate the approx i mate location of the spiral arms detected 
by IFukagawa et alJ d2004l) . The white curve is a possible in- 



ward spiral arm revealed by this observation. The outer ring, 
at 350 AU, indicates the measured outer radius of the contin- 
uum emission. The inner ring, at 280 AU, is the location of the 
feature observed by Pant in et al J (12005) at 20.5 jjm. 




"0 




Fig. 5. Images of the continuum 2.8 and 1.4 mm emission of 
AB Aur at 1.4 mm. 1.4 mm image: resolution 1.3", contour 
spacing 1.3 mJy/beam (22 mK, 2 cr) 2.8 mm image: resolution 
2.0", contour spacing 0.6 mJy/beam (16 mK, 2 cr) The ellipse, 
at 350 AU, indicates the measured outer radius of the contin- 
uum emission. 



radius is found at 1.4 mm than in the CO lines (see Table 1). 
This is presumably due to a combination of the smaller dust 
opacity with a progressive fall of the surface density inside the 
ring. When modelled with a sharp edge, this naturally results in 
the most optically thin tracer, namely the dust emission, having 
the largest apparent inner radius. 

A simultaneous fit of the 2.7 and 1.4 mm continuum 
data allows to derive the spectral index of the dust emis- 
sivity, /3 = 1.4 ± 0.2, assuming the emissivity k(v) = 
O.lcm^g-V/lO^Hzy 3 . 



6 



Pietu, Guilloteau & Dutrey: Sub-arcsecond images of AB Auriga 



3.3. Line Emission: the Disk 

The double peak profiles, combined with the strong velocity 
gradient along the major axis of the emission which is found 
in all transitions, suggest that the material around AB Aur is 
dominated by a circumstellar disk in Keplerian rotation. 

We thus attempted to model the emission by a sim- 
ple Keplerian disk model, in hydrostatic equilibrium, as 
perf ormed successfu l ly fo r T Tauri s tars like DM Tau 
ilGuilloteau & Dutrevl [jUj) , GM Aur dDutrev et all Il998l) 
and other PMS stars including the Herbig Ae star MWC 
480 ( Si mon et alJ l2000). We modelled the circumstellar mat- 
ter, gas and dust, usi n g the same strategy than those ap- 
plied in iDartois et all J2003b for the DM Tau disk. The 
12 CO J=2— > 1 transition traces the CO surface of the disk while 
the 13 CO J=l— >0 & J=2— >1 lines trace material closer to the 
mid-plane. 

Clearly, AB Aur is surrounded by a very large disk. All 4 
spectral lines give a consistent value for the orientation of the 
disk axis, -30° ±1° East from North. The outer radius is of 
order 1000 AU, except for C ls O J=l— >0 which is significantly 
smaller (600 AU). Note that the J=2— >1 transitions, specially 
from 12 CO may underestimate the true radius because of a lack 
of sufficient short spacings. 

Both 13 CO transitions give a consistent value for the rota- 
tion velocity at 100 AU, vsin(0 = 1.60 ± 0.05 km.s -1 , and a 
slightly different value from 12 CO. When attempting to fit the 
data with a simple Keplerian disk model, the derived inclina- 
tions are different (34, 36 and 43°), and the derived stellar mass 
(0.9 -1.2 M Q ) far from the expected value for an AO star. From 
Table 1, it is clear that a simple Keplerian disk model fails to 
fit properly all observed transitions, even if account is taken of 
a possible vertical temper ature gradient in the disk as has been 
observed in DM Tau by IDartois et alJ J2003I) . The most sig- 
nificant result is that the velocity law deviates from Keplerian 
motions. From 13 CO J=2— >1, the exponent is v = 0.42 + 0.01, 
instead of v = 0.50 for Keplerian motions, and the J=l— >0 
gives a consistent result. The law deduced from I2 CO is much 
steeper, with an exponent of order 0.8. 

It is of course conceivable that the 12 COJ=2— >1 data be 
discarded in such an analysis, because of the confusion with 
the surrounding molecular cloud. Indeed, the fit by emission 
from a disk leaves significant residual, and cannot represent for 
example the excess emission in the N-E which is seen in CO 
at velocities between 4.7 and 6.8 km.s~'. The large turbulent 
width found in this analysis is just a result of attempting to 
fit non-symmetric features with a pure rotating disk. Despite 
these caveats, an important result from the CO data is the high 
temperature (Tioo - 65 - 70 K) required in the inner 100 AU. 

The two transitions of the 13 CO isotopologues yield rea- 
sonably consistent results. Some of the differences can be at- 
tributed to the importance of the emission coming from the 
continuum spiral in the 13 COJ=2— >1 data. This bright emis- 
sion dominates the total flux, and biases the derived inclina- 
tion towards higher value, because the spiral curves inwards 
towards the North- West. If the aspect ratio of the I3 CO J=2->1 
integrated emissio n is taken as a measure of the inclination of a 
simple disk, as did Ma nnings & Sargenll I I 19971) . quite high val- 



ues are found (^ 70°). Significantly lower values are obtained 
when the constraints from the kinematics of a rotating disk is 
incorporated, but the result on the inclination is still biased be- 
cause of the emission from the spiral. The 13 COJ=l— >0 line 
being observed with lower angular resolution is less sensitive 
to the spiral feature, and yields an inclination i = 39 + 2°. Given 
the departure from circular symmetry, lower values (e.g. down 
to 25°) cannot be excluded. We note also that it is necessary to 
introduce an inner radius of about 70 AU to best reproduce the 
I3 CO emission (see Table). This inner gap is clearly visible in 
the integrated intensity maps of FigQ 

The most clear result from the 13 CO analysis is that the 
best fit to the velocity law is obtained with an exponent v = 
0.41 +0.01. This is a very robust result: we checked that the de- 
rived exponent is independent of the assumed inclination over 
a range from 20° up to 50°. This robustness also applies to 
most other parameters (column density law, inner and outer 
radii, temperature law, turbulent width). Using a stellar mass 
M, = 2.2 M Q for AB Aur, the Keplerian speed at 100 AU 
is =i 4.4km.s~ 1 . Accordingly, for the nominal inclination of 
39°, the orbital velocities are sub-Keplerian (v(r) < y/GM t /r) 
out to the disk edge. However, if the inclination is lower than 

22°, the orbital velocities are super-Keplerian outwards of 

100 - 150 AU. 

C 18 J=l— >0 data, although more noisy, also give a i 39° 
inclination. This result is important, since C ls O J=l— >0 is the 
most optically thin of all observed transitions, and thus the less 
likely to be affected by an extended envelope. The velocity 
exponent found from C ls O v = 0.47 ± 0.04 (which still re- 
mains compatible with the 13 CO results at 2<x), is essentially 
Keplerian, but the projected velocity at 100 AU is much smaller 
than expected from an AO star, as in the case of other transi- 
tions. To reconcile the velocities with the expected stellar mass 
would require an inclination of ^ 23°. Note also that, because 
of its optical thinness, this transition mostly samples the inner 
radii with little contribution from the outer parts of the disk 
(r > 300 AU) and the envelope. Also, an independent con- 
straint on the surface density and the temperature laws is no 
longer possible; for an optically thin J=l-0 line, the data only 
const rain the ratio E/T and the value of p - q ( Da rtois et alJ 
120031) . 

For all transitions, we have also checked whether adding a 
radial velocity component would provide a better fit to the kine- 
matic pattern. For 13 CO and C 18 , the upper limit on a radial 
velocity is below 0.05 km.s . For 12 CO , a slightly better fit 
is found by adding outward motions of order 0.2 km.s 1 . The 
significance of this result should not be overestimated, since 
the excess emission towards the North near Vlsr — 6.5 km.s -1 
(see Fig|2j may bias the result. Note that we can discriminate 
between inward and outward motions only because we know 
the full 3-D orientation of the AB Aur disk (the Sou thern part 
being towards us, in agreement with Fu kagawa et alJ (12004)). 

4. Discussion 

The observations presented above bring new information about 
the circumstellar environment of AB Aur, which appears to be 
very different from other proto-planetary disks. 
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Table 1. Best parameters for the AB Aurigae disk derived from^ 2 minimization in the UV plane. 



Assumed Distance 


D(pc) = 140 


Lines 


12 COJ=2^1 t 


13 COJ=2^1 X 


13 COJ=1^0 $ 


C 18 OJ=1^0 ft 


Systemic velocity 
Orientation 
Inclination 
Inner radius 
Outer radius 
Turbulent linewidth 


V LSR (km.s- 1 )= 5.89 ±0.01 
PA O = -27.5 ±0.5 
i(°)= 33 ±1 
R m (AU) = 45 ±3 
R out (AU) = 1050 ± 10 
Av(km.s-')= 0.38 ±0.02 


5.84 ±0.01 
-31.3 ±0.3 
42 ±1 
72 ±5 
890 ± 10 
0.22 ±0.02 


5.87 ± 0.03 

-30 ±1 

39 ±2 

77 ±5 
1300 ±100 
0.26 ±0.02 


5.84 ±0.02 
-29 ±1 
36 ±3 
67 ±6 
600 ±60 
0.18 ±0.05 


Column Density 

j_ i r\ r\ att 

at 100 AU 
exponent 


Molecular Column d 

E 100 (cm- 2 ) = (1.4 10 19 ) 
P = (2.5) 


snsity law: Z(r) = Xiooi 

2.3 10 17 ±0.3 10 17 
2.4 ±0.1 


r y-t> 

100AU ' 

3.2 10 17 ±0.8 10 17 
2.5 ±0.3 


9.5 10 16 ±1.3 10 16 


Temperature 
at 100 AU 
exponent 


Temperature la 

r 100 (K) =68 ±1 

q = 0.77 ±0.02 


w: T(r) = T m ( 100 r AU 

38 ±2 
0.16 ±0.03 


34 ±2 
0.10 ±0.14 


20 ±3 
-0.6 ±0.2 


Velocity at 100 AU 
Velocity exponent 


Velocity law: 
V 100 (km.s- 1 )= 3.06 ±0.06 
v = 0.82 ±0.02 


VW = V 10 o(to^u)- 
2.54 ±0.03 
0.42 ±0.01 


V 

2.73 ±0.12 
0.37 ±0.02 


2.75 ±0.13 
0.47 ±0.04 



$ ) For the l3 CO lines, the error bars are derived from the l<x formal 
errors from the x 1 fit- We checked that all parameters (except of course 
V*ioo) are essentially independent of the inclinations over the 25 - 45° 
range (see Section l33l . 

t) For I2 CO, multiple minima exist. We have selected a solution in 
which the scale height is determined from the 13 CO results: /jjoo = 12 AU, 
h = -1.40. The error bars, when mentioned, represent the typical (1 
cr) variations of the parameters obtained when varying all other ones. 
Numbers in parenthesis represent assumed (fixed) values. 

ft) For the C 18 0, only p - q is actually constrained: p — q = 3.6 ± 0.2. 
See text for details. 

* ) For the continuum data, the temperature was taken as 35 K at 100 AU 
with an exponent of q = 0.3. The surface densities and mass (and their 
errors) are given for p = 1.4. An additional 30 % uncertainty should be 
added for these quantities, because of the error on [}. 



Continuum Results * 
Dust: k v = k„ x (t^tt^ 


f 




Absorption law 


k (cm 2 /g) = 




(0.1) 


Dust exponent 


P = 




1.4 ±0.2 


Surface Density law: £(r) = Xioo( 


100 AU ' 


Surface Density 


Eioo (H 2 cnr 2 ) = 


6.3 10 23 ±0.8 10 23 


at 100 AU 


S 100 (g.cirT 2 ) = 




2.7 ±0.3 


exponent 


P = 




2.3 ±0.2 


Inner radius 


Ri (AU) = 




115 ±5 


Outer radius 


Rj (AU) = 




350 ±30 


Mass 


M rf (M ) = 


1.6 10~ 2 ±0.2 10~ 2 



4.1. The "spiral" 

First, the sub-arcsecond image of the continuum emission at 
1 .4 mm extend s inwards the s piral structure detected in the near 
IR by Fukaga wa et alJ J2004I) . Considering the sense of rota- 
tion, the new portion of the spiral is trailing, as the previous 
arms. Dust emission at long wavelengths being proportional to 
KyLT, this enhanced brightness could be due to 3 effects: 1) 
a column density, 2) a change in dust properties, with larger 
grains in the spiral, or 3) a larger temperature. However, grain 
growth cannot explain the enhanced molecular emission which 
is also seen in the spiral. It is not straightforward to disentan- 
gle between the two remaining effects from the present data 
only. Taken at face value, the analysis of 13 CO lines seem to 
indicate a flat (or even rising) temperature throughout the disk. 



However, this result is to be taken with care, since the 1J CO 
lines are optically thick only in the inner spiral, whose struc- 
ture is not adequately fitted by the disk model. A better in- 
sight is given by the fa ct that the C 18 OJ- l— >0 line emission 
is proportional to £/T dDartois et al. 2003), rather than YT for 
optically thin dust emission. Since the brightness contrast in 
C 18 OJ=l— »0 is similar to that in the 1.4 mm continuum, it 
thus suggests the enhanced brightness is a column density ef- 
fect. 

The difference between the apparent inner radius in dust 
(110 AU) and in CO (- 70 AU from all lines) indicates 
that the structure of the inner region is much more complex 
than a simple hole, in sharp contrast with t he case of GG 
Tau JPutrev et alJl994lGuilloteau et alJl999l) . but its detailed 
structure is beyond the resolving power of these observations. 
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4.2. The Inclination Problem and non-Keplerian 
Velocities 

By revealing non circular structures, these observations offer 
some explanation for the widely different inclinations which 
were derived in the past from mm tracers and optical images. 
High inclinations are ruled out by our data. Although an incli- 
nation of =i 38° seems adequate to represent most of the emis- 
sion, it is certainly biased by the non-circularity of the emis- 
sion. Given the uncertainties, inclinations as low as 23° cannot 
be excluded. 

On the other end, we note that the inclinations derived for 
the inner disk (< 2 AU) through near-IR interferometry by 
lEisner et alJ J2003I) could be biased towards low values by non- 
circular structure, if the "spiral" pattern persists at such dis- 
tances from the star. 

In any case, the orbital velocities significantly differ from 
the Keplerian speed. In addition, the turbulence (see Table 1) 
appears somewhat larger in AB Aur (0.2 km.s~' if we exclude 
12 CO) than in other sources like e.g. DM Tau. Unfortunately, 
the exact value of the inclination plays a major role in the un- 
derstanding of the environment, because depending on the as- 
sumption about projection effects, the orbital velocities of the 
material surrounding AB Aur can be sub-Keplerian or super- 
Keplerian throughout most of the disk. 

4.3. Kinetic Temperature 

In this respect, the AB Auriga disk behaves like those around 
other T Tauri stars. There is a clear evidence for a kinetic tem- 
perature gradient from the disk mid-plane to the disk surface 
and envelope, with the disk plane being cooler (around 30 K, 
as traced by 13 CO and C lo O) than the envelope (~ 70 K at 100 
AU, as traced by 12 CO, see Table 0, as expected from heat- 
ing by the central star. The temperature remains high (^ 30 K) 
throughout the disk, as expected from the much higher lumi- 
nosity of AB Aur compared to typical T Tauri stars. However, 
the temperature derived from I2 CO falls below that derived 
from 13 CO lines beyond a radius about 400 - 500 AU. This 
could be due to self-absorption by the (colder) envelope. 

4.4. Molecular Abundances and Disk Mass 

The derived outer radius of the dust disk is much smaller 
(=! 400 AU) than that the molecular disk (^ 1300 AU). At 200 
AU, the ratios of the molecular column densities and H2 sur- 
face density derived from the dust indicate a molecular abun- 
dance [ 13 CO/ H?1 ^ 5 10~ 7 , similar to that of the Taurus cloud 
nO- 6 . lFrerking et alJ JlQ82h . The value of the dust opacity in- 
dex j3 = 1.4 + 0.2 (giving K220GHZ = 0.012 cm 2 /g), is at the 
high end of the values found in circumstellar disks (and in par- 
ticular of those observed simultaneously with A B Aur, w hich 
share the same flux calibration, see Section 2.1). lDutrev et alJ 
(1996) indicate a mean value of f3 — 0.95 + 0.11 from an en- 
semble of disks. This suggests that the disk around AB Aur 
may be younger than average, and that dust grains have not 
evolved as much as in disks around T Tauri stars. The rela- 
tively high temperatures (above 30 K) prevent condensation 



on dust grains, and probably explain why CO is not signifi- 
cantly depleted, con trary to most other circumstellar disks (e.g. 
lDutrevetalJdl996l) 1. 

The lack of detectable dust emission beyond 350 AU could 
indicate that the dust opacity drops there. Since the molecular 
content only changes smoothly with radius, this suggest that 
the dust absorption coefficient at 1 .4 mm is smaller beyond 350 
AU than inside, i.e. that the dust is significantly less evolved 
outside. Another alternative is a steepening of the surface den- 
sity distribution beyond 300 AU, since the values derived 
from the dust (p = 2.3) is slightly lower than that derived from 
the CO molecules (p - 2.7). The latter possibility is however 
not supported by the C 18 data, which are only sensitive to 
the inner regions. It is tempting to relate the change of dust 
prop erties near 350 AU to the bright 20 fim ring detected by 
iBouwman & Pantinl J2003I) at 280 AU. 

Both the high f3 value, and the possible dust opacity drop 
beyond 350 AU suggest that the dust in AB Aur is less evolved 
than in other T Tauri disks. 

5. Origin of the disturbances 

The AB Aur environment, although dominated by a rotat- 
ing disk, is clearly highly disturbed and far from the quasi- 
equilibrium Keplerian stage encountered around other young 
stars. We discuss here some mechanisms which could lead to 
these peculiarities. 

5.1. Multiplicity ? 

Ring-like structures are easily formed by tidal trun cations in bi- 
nary systems, as illustrated b y the case of GG Tau JDutrev et alJ 
ll994UGuilloteau et al. 1999). AB Aur has never been reported 
as having a companion. Limits on any companion mass can 
be derived from the liter ature. For example, us ing the bright- 
ness profile published by (Fukaaawa et al. 2004) in t he H b and, 
and the evolution models from Baraf fe et alJ 41998, 2002), we 
find a (conservative, especially at large radii) upper limit on 
any coeval companion of less than 0.02 M in the field of view 
of the Subaru image. This upper limit on the mass holds for 
120 < r < 1500 AU. For closer objec ts, the limits are more 
difficult to quantify Jpirzkal et alJ dl997h give an upper limit of 
0.25 Mq down to . 4" (6 AU). A stronger limit can be derived 
from lLeinertetalJ ( 1 1994 in the range 0.07 - 1.0" (10 to 140 
AU), as their speckle observations show that the flux of any 
companion in K band could not exceed 1 to 3 % of the flux of 
AB Aur, for this range of separations. This gives an upper limit 
of about 5.2 - 6 to the K magnitude of a possible companion, 
which translates into a mass limit of 0.05 - 0.3 M for ages 
between 1 and 10 Myr. Unfortunately, a companion of mass 
M c <K M„ on an orbit with semi-major axis a could create a 
gap of half-width w a few times the Hill's radius 

where F is an empirical factor of order 2 - 3(Morbidelli 2002). 
If viscosity is not negligible, the conditions for gap creation 
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and the resulting gap width are diffe rent. U sing an a prescrip- 
tion for the viscositv. lTakeuchi et alJ i 19961) indicate a gap half 
width 

w^l.3aA 1/3 (2) 
where A is the ratio of strength of tidal effects to viscous effects 



\M») 3a(h/r) 2 



(3) 



Using h/r 0.1 as derived from our measurements, and 
a = 0.01, we find A > 1, for M c > 0.05 M . Hence, viscos- 
ity is unlikely to affect our conclusions. Accordingly, even a 
0.05 M Q companion orbiting AB Aur around 40 AU could con- 
ceivably evacuate the inner hole. Whether it could sustain the 
non Keplerian velocities requires detailed studies. Higher sen- 
sitivity optical observations are required to setup more stringent 
limits on the existence of any companion. 

5.2. A gravitationally unstable disk ? 

Table^shows that the total mass of the system made by the disk 
and the envelope is not important enough to be self-gravitating. 
This can be better assessed using Toomre's Q criterium: 



Q = 



ttGE 



< 1 



(4) 



where c s is the sound speed, and Q the angular momentum. 
With our power laws, Q is also a power law of the radius Q(r) = 
Qoir/Ro)-' with 



Go = — / A ^" V ° 

t = 1 + q/2 + v - p 



(5) 
(6) 



Using the values from Tabled we find go 11 at 100 AU 
and t = -0.7 + 0.2, well above the stability limit. To have 
the disk massive enough for instability would only be possi- 
ble if the CO depletion was high (> 20) and the dust opacity 
10 times lower than we assumed. Both possibilities are highly 
unlikely as discussed in Section l4~4l For example, lowering the 
dust opacity by a factor 10 would require a value of = 3.9, as- 
suming the absorption coefficient at 10 12 Hz remains constant. 
Furthermore, t being negative, the stability increases with dis- 
tance from the star. 

5.3. Possible encounters 

Another possible explanation is a past encounter with a star. 
We have been looking in the SIMBAD database for all stars in 
a radius of 30 arc-minutes centered around AB Aur. Of a total 
of 31 stars referenced, 16 of them possess prop er motion mea - 
surements. 9 come from the survey of Jones & Herbig ( 1979). 
4 com e from the Tycho Reference Catalog ilRoeser & Bastian 
1988), and 3 Hipparcos measurements were av ailable (4 all to- 
gether with AB Aur) - JPerrvman et aljfl~997h . A very simple 
analysis shows that at least 2 of them could have encountered 
AB Aurigae in the past. JH433 could have encountered AB Aur 
at any time older than 35 000 years, whereas there could have 



been an encounter between AB Aur and RW Aur, some 500 
000 years ago. We should also mention that according to the 
early measurement of Jones and Herbig, SU Aur could have 
met AB Aur is a close past, but that Hipparcos measurements 
rule out this possibility. 

5.4. A young disk ? 

Rather than being affected by multiplicity, the situation may 
suggest instead that AB Aur is surrounded by a young circum- 
stellar disk which has not yet relaxed to the Keplerian stage. 

In the formation of a star + disk from the collapse of a ro- 
tating interstellar cloud, the low specific angular momentum 
directly accretes to form the star, but the higher angular mo- 
mentum accretes through a disk. ICassen & Moosmanl dl 98 lit 
have pointed out that, because of the projection of the mo- 
mentum along the disk axis, the accreting material always have 
lower angular momentum than the Keplerian value where the 
accreting material trajectories cross the disk. Ho wever, when a 
comp lete history of the accretion is incorporated. lStahler et alJ 
(1994) have shown that the disk evolution, which occurs be- 
cause of the drag force induced by the accreting material, re- 
sults in a semi-steady state situation in which the outer disk 
is super Keplerian. Bot h studies concl ude that a massive ring 
should be formed. Stahle r et al.l ll 19941) show that this ring ac- 
cumulates mass, and acts as a temporary storage for the excess 
angular momentum brought by the accreting gas. 

It is tempting to identify this ring with the enhanced density 
region which is seen in continuum and line at about 150 AU 
from the star. If the inclination of the disk is low, 23°, the 
measured rotation velocities exceed the Keplerian speed at all 
distances above 100 AU from the star (our measurements are 
insensitive to the exponent of the rotation velocity within the 
inner 100 AU). The apparent spiral arms could be the result of 
instabilities which necessarily occur in such a situation, since 
the ring is a transient phenomenon. A possible difficulty in this 
interpretation resides in the timescales over which the ring and 
supra-Keplerian velocities persist. Accretion on the disk only 
occurs over the infall timescale, which remains to first order 
similar to the free-fall time ts, i.e. a few 10 4 years at most. 
The age of AB Aur, f AR Aur is estimated t o ran ge between 1 
and 4 Myr. 4 Myr is cited bv lDeWarf et alJ J2003I) on the basis 
of a coevality with the nearby star SU Aur, which has similar 
high proper motions to AB Aur. Lower ages are usually quoted 
by other authors, AB Aur being among the youngest known 
HAeBe star. However, AB Aur is already relatively hot, being 
of spectral type AO V, so that an extremely young age, a few 10 5 
years, appears to be ruled out from the currently available evo- 
lutionary tracks. The main accretion phase on the disk should 
thus likely be finished. 

However, the longest timescale over which the circumstel- 
lar disk evolves is the viscous timescale, 



fvisO) 



r 

v Xmb 



ac s H(r) 

in the a prescription of the turbulence. We find 



akT(r) 



fvi s (100 AUXr/lOOAU) 1 



(7) 



(8) 
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Using a = 0.01, and the temperature and velocity from Table 
1, we find ? v ; s 10 6 yr at 100 AU, which is comparable to the 
estimated age of AB Aur. It is thus conceivable that the AB 
Aur disk has not relaxed to a purely Keplerian stage and stills 
exhibit a disturbed kinematic and density pattern. However, it 
is important to mention that the viscous timescale strongly de- 
pends on the assumed viscosity parameter a, but little on the 
stellar properties. As a consequence, the above value is also 
applicable to T Tauri disks which exhibit clear Keplerian signa- 
tures, like those around DM Tau, LkCa 15 or GM Aur. Unless 
the viscosity is different in both type of objects, this suggests 
that AB Aur is younger. 

6. Summary 

We have shown from observations of dust and 4 lines from the 
CO isotopologues that AB Aur is surrounded by a flattened, 
non-circularly symmetric disk, in rotation around the star. As 
expected from the high stellar luminosity, the disk is warm, and 
CO is not significantly depleted. The determination of the disk 
inclination is biased by the non-circularity of the brightness 
distribution, but the rotation has clear non Keplerian character- 
istics. No completely satisfactory explanation exists for these 
unusual characteristics. The disk does not appear to be massive 
enough to be self-gravitating. Upper limits on any companion 
are not stringent enough to conclude. No evidence for infall 
motions could be detected. It is tempting to see AB Aur as an 
example of a very young object in which traces of the accre- 
tion from a rotating envelope onto the disk surface have not yet 
been erased by the dynamical and viscous evolution of the disk. 
An independent indication of youth comes from the dust prop- 
erties, since the /? index may indicate that dust grains have not 
evolved in the AB Aur disk as much as in other proto-planetary 
disks. 
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